Dispersive Manipulation of Paired Superconducting Qubits 
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We combine the ideas of qubit encoding and dispersive dynamics to enable robust and easy quantum infor- 
mation processing (QIP) on paired superconducting charge boxes sharing a common bias lead. We establish a 
decoherence free subspace on these and introduce universal gates by dispersive interaction with a LC resonator 
and inductive couplings between the encoded qubits. These gates preserve the code space and only require the 
established local symmetry and the control of the voltage bias. 
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Superconducting nano-circuits consisting of charge boxes 
(CB) 01 are among the most promising candidates for a quan- 
tum computer. Coherent control of a single charge qubit 0], 
long decoherence time 0] and, more recently, coupled two 
qubit systems 0] have been demonstrated. But despite this 
encouraging experimental progress, there are serious difficul- 
ties with superconducting QIP which may appear insurmount- 
able. The first is the severe decoherence experienced by these 
macroscopic qubits, which are coupled to a large number of 
degrees of freedom in their environment including control cir- 
cuitry jsl . A few methods have been employed to reduce the 
decoherence 00], but they usually require sophisticated ma- 
nipulation or significant overhead in the control circuitry. The 
second major difficulty comes from the imperfect control re- 
alizable in solid state systems. Specifically, one finds it diffi- 
cult to achieve controllable couplings between superconduct- 
ing qubits, since the commonly used hard-wired inductive or 
capacitive couplings are untunable. Great effort has been exer- 
cised to realize controllable couplings |8, 9, 10]. Schemes al- 
lowing to compute with invariable couplings were also studied 
fTll[l2il . Others have recently discussed to use an LC circuit 
to actively mediate the interaction between superconducting 
qubits H3til4l1 . 

The requirement to reduce decoherence and the desire for 
the easiest manipulation apply to all QIP implementations. 
Unfortunately, it is not always easy to accomplish both - ac- 
tually the goals are often contradictory - since reducing deco- 
herence may require extra complication in the manipulation. 

In this work we show how to achieve both goals. Using 
a closely placed pair of charge boxes (PCB) sharing a com- 
mon bias lead as the logic qubit, we can encode informa- 
tion in a fashion immune to collective noise, which is the 
dominating decoherence source in our setting. We introduce 
LC resonators inductively coupled to the PCBs whose vir- 
tual excitations allow us to manipulate the PCB dispersively; 
all interactions will be off resonance, without energy trans- 
fer Olil, and thus a logical qubit stays within its encoding 
space even during manipulation. By inductively coupling the 



CBs and taking advantage of dispersive dynamics again, con- 
trolled phases can be induced between logical qubits. 

Combining dispersive dynamics and encoding offers a new 
method for QIP. It overcomes the major difficulties of super- 
conducting CBs in a realistic and very simple fashion. The 
only control required after initialization is to change the volt- 
age bias of the PCB. Another advantage of our method is that 
it relies only on noise symmetry over short distance and so 
this is a realistic technique for scalable QIP over large sys- 
tems. Though we focus on superconducting QIP in this paper, 
the general principle of dispersive manipulation of encoded 
qubits will be valuable for many other QIP systems. 

Charge qubits and the dominating noise. - In its simplest 
form, the charge qubit is just a superconducting island volt- 
age biased through a Josephson junction. The Hamiltonian 
for the CB system is Ec{n — rig)"^ — Ej cos ip Ol, where the 
charging energy Ec = (2e)-^/2Ct, Ct being the total capaci- 
tance of the island, is much greater than the Josephson energy 
Ej ~ Ic^o/2tt, and (p is the phase drop across the junction. 
When biased close to Ug = CgVg/2e =1/2 (Cg is the gate 
capacitance), it provides an effective two-state system which 
can be used as a qubit. In the spin 1/2 notation. 
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where the spin up or down states correspond to n=0 or n=l 
excess Cooper-pairs on the CB. The effective field i?^ ~ 
Ec{2ng — 1) can be tuned by changing the gate voltage Vg. 
In order to control = Ej, we use a flux biased small dc- 
SQUID to replace the junction, whose critical current is max- 
imal {Ic = Ic) when the flux bias is off and vanishes (/c = 0) 
when it is $o/2 = 

The dominating decoherence sources in this system are cir- 
cuit noise in the voltage bias 01 and charge fluctuations in the 
background (known as "charge noise") jlTlllSll . as indicated 
in Fig. ^(a). The circuit noise is described by the well known 
"Spin-Boson" model The charge noise is less well under- 
stood, but it is now generally believed to be caused by fluc- 
tuations of the impurity charges in the substrate flTll . Taking 
into account the noise sources, we have the total Hamiltonian 
for the system as if = Hcb + Hz + Hb + Hint, where 
Hz and Hb are the Hamiltonian of the circuit fluctuations. 
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modelled as a collection of harmonic oscillators ^T^, and the 
Hamiltonian of the background charge. Since these noise per- 
turb the voltage bias of the CB, the interaction Hamiltonian 
Hint has the form Hint = cr^E^ + a^E^, where E^ and 
E^ are environment operators (with the coupling strengths 
included). Here we focus on the fluctuations in the voltage 
bias, the dominating source of decoherence and neglect noise 
in the field (the critical current), as practised customarily 
fisl fl9il . This is because charge qubits are insensitive to 
flux noise and the effect of the fluctuation in the suppres- 
sion field is secondary to the bias voltage variations discussed 
above. Detailed treatment of the two voltage noise sources 
and their influence on the charge qubit system can be found in 
the literature fl^fioll . For our purpose, the nature of the en- 
vironment and specific form of E^ and E^ are not essential, 
therefore we do not go into detail here. 
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FIG. 1: (a) A pair of closely spaced CBs sharing a common lead 
used as the encoded qubit. The SQUIDS of both CBs are induc- 
tively coupled to an LC resonator whose virtual excitation allows the 
two charge boxes to exchange energy quanta. The circuit noise and 
background charge fluctuations are represented by Z and B respec- 
tively, (b) An inductively coupled PCB array. For clarity the PCBs 
are schematically represented by two boxes, the inductive coupling 
is understood to be between the dc-SQUIDS of the PCBs. 

Paired charge boxes and DPS encoding. - As shown in Fig. 
[l](a), we use two capacitively coupled identical CBs (a and b) 
with a common bias-lead as an encoded qubit. The small ca- 
pacitive coupling, Cc Ct, is not essential for the encoding, 
but necessary for the encoded two qubit gates. The Hamil- 
tonian of the PCB system is HpcB = Yl,i=a b l^^^'^t ^ 
^i'^i) ^ l^l^b' where 7 = ^r^c is the coupling energy. 

Since the two CBs share the same lead, obviously they are 
biased at the same voltage and they will experience the same 
circuit noise. In addition the nano-scale charge islands are put 
very close to each other Therefore, they will experience the 
same charge fluctuations too (more discussion on this point 
will be given later). Hence the CBs experience "collective 
decoherence," meaning the noise sources couple symmetri- 
cally to them, which naturally gives rise to "decoherence- 
free encoding." For a review of decoherence-free subspace, 
see 1I20I1 and references therein. Here we have the simplest 
case of the DFS, with |0) = | ia]b) and |1) = | Taib) as 
the decoherence-free logical states. The way this works can 



easily be seen: as a consequence of the collective decoher- 
ence the coupling between the PCB and the environment is 
(cr^ + a^) * [E^ + E^), which annihilates the two logical 
states given above. Therefore the PCB system will not get 
entangled with the environment if it is initialized and kept in 
the DFS. Physically, the encoded qubit' s immunity to noise 
stems from the fact that the CBs acquire random but opposite 
phases. 

To prepare the system in the DFS, we bias the PCB far off 
the degeneracy point Ug = 1/2. In the spin- 1/2 picture, this 
corresponds to applying a strong field in the z direction. At 
low temperatures, the spins will line up with the field, and the 
PCB relaxes to the state | [aib)- Keeping both off, we 
change (same for a and b) to —27. This cancels the bias of 
aoTib making its total B^^^ = _B- +27 = 0. After that we turn 
on some B^. After a time n/B^, the state of b will become 

tb), and the system is prepared in the DFS state | jaTb)- 
Now we turn off Bf as well, and from now on the B^ fields 
for both CBs will remain off, by biasing the dc-SQUIDS of the 
PCB at $0/2. Since B^ fields will remain off and need not 
be tuned after the initialization, the leads tuning Ic of the PCB 
could be heavily filtered to keep out the noise once the system 
is initialized. Alternatively we could make use of the noise 
free constant flux-bias techniques such as that demonstrated 
in 12 111 . In practice it can be difficult to suppress Ic of the 
PCB precisely to due to the finite self inductance of the dc- 
SQUID. However as shown in f25| if low self inductance dc- 
SQUIDS are used (^ «: 1) the B^ field at $o/2 bias point is 
negligibly small compared to B~ fields used for computation 
and can thus be safely dropped. Many schemes [ll |^ |21 rely 
on this fact too. One notices that the logical states |0) and 
|1) are always degenerate regardless of the voltage bias Ug, 
therefore there is no evolution in the idle mode, regardless of 
the voltage bias or noise in it. To readout the state of the PCB, 
a measurement of its CB a or 6 will suffice, which is readily 
accomplished with developed techniques 1 1 ') | . 

As seen above, in realizing DFS encoding with the PCB, 
we lose considerable freedom in manipulating the system. 
First, in order to guarantee symmetrical coupling to the cir- 
cuit noise, the two CBs share the same lead and hence they 
are always biased at the same voltage. More importantly, we 
must ensure that operations on the PCB do not drive the sys- 
tem out of the DFS, otherwise the immunity to noise is lost 
I20I1 . This is why we must keep the B^ fields for the CBs off: 
they flip the states of a single CB and hence do not preserve 
the DFS. Therefore the only control left is the voltage bias of 
the PCB, which clearly is insufficient for universal QIP on the 
PCB. To deal with this difficulty, we introduce a LC resonator 
inductively coupled to the PCB system: 

The LC resonator inductively coupled to the PCB. - As 
shown in Fig.Q](a), we inductively couple the dc-SQUIDS of 
the two CBs in the PCB system symmetrically to an initially 
unexcited LC circuit. Even in the ground state its vacuum 
fluctuations bias the dc-SQUIDS of the PCB off $o/2 making 
charge tunnelling possible. The Hamiltonian for the PCB-LC 
system is H = HpcB + Hlc 



Hcoup, where Hlc = uja^a. 
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and L refer to frequency and inductance of the LC-resonator; 
M is the mutual inductance between SQUID and resonator. 
When the PCB and LC are far off resonance, the effect of 
the LC can be neglected. On the other hand, when we tune 
the bias of the PCB such that it is close to being in reso- 
nance with the LC resonator, within the framework of Ro- 
tating Wave approximation the above Hamiltonian becomes 
H = HpcB + Hlc - i9Y.i=a,bi^t o- - crra^), where 
CT^ = (cT^ ± ia^)/2 are the ladder operators. Notice that 
Hpc B contains a coupling term , in contrast to the standard 
Jaynes Cummings Model which modifies the dynamics sig- 
nificantly, as will be seen below. 

Let S = — ujhe the detuning. If we let the PCB and the 
(initially un-excited) LC resonator interact right in resonance, 
i.e., |J - 27I < 5 < w or |(5 + 27I < 5 < w Q, state 
transfer occurs between the PCB and the LC resonator fl3l 
Hill . This is not allowed in our scheme, since it will drive 
the PCB out of the DFS. Besides, once the LC resonator is 
excited, we will have additional decoherence due to the finite 
quality of the LC resonator iHslH . Therefore, we only work 
in the dispersive region, g ^ \5 zt 2j\ ^ uj. In this case, the 
PCB and the LC resonator cannot exchange energy because 
of the large detuning. However the virtual excitation of the 
LC resonator gives rise to an effective interaction between the 
two CBs iH, 



b) in the PCB. Then when we tune the voltage bias of the 
PCB such that it interacts with this LC resonator dispersively 
((7 <C 1^' ± 27I <C i^'), a phase gate will be obtained due to 



H, 



where the first term describes the Stark shift, which can be 
neglected since it is the same for both logical states, and the 
second term is the effective exchange interaction caused by 
the exchange of a virtual photon. It preserves the DFS (it 
changes | lalb) to | iaib) and vice versa) and acts as a log- 
ical X gate on the PCB. Starting from |0) = | iatb). let- 
ting the PCB evolve under the effective Hamiltonian (|2ji for 
a time t = '^^^'^^^'''> or |, we can swap the states of a and b 
or generate a maximally entangled state between them. The 
LC resonator is initially in the vacuum state and will not be 
excited due to the dispersive interaction with the PCB; there- 
fore unlike in previous schemes flslO we are not subject to 
decoherence caused by its finite quality. Also, notice that our 
dispersive scheme is fundamentally different from the previ- 
ous method of using the virtual excitation of a large LC circuit 
capacitively coupled to all the qubits fsl, in which the LC- 
frequency is much larger than the CB-energies and the cou- 
pling strength is tuned by changing the Ej's of the qubits. In 
our scheme the two energies are close (though the detuning 
is large) and the -Ej's are always off; the coupling strength 
is controlled by simply changing the detuning. Our scheme 
takes advantage of the quantum exchange effect (Eq. (2)) as- 
sisted bya virtually excited quantum state of the LC-qubit 
system I24l IZSll . It offers noise protection with simplest oper- 
ation using very realistic parameters (see below), which was 
not available in previous schemes. 

To realize SU{2) on the PCB, we still need a phase gate. 
This can be accomplished by using another LC resonator, not 
depicted in Fig. ^{a) and at a frequency ui' very different 
from uj, inductively coupled to the dc-SQUID of only a (or 
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the Stark shift 



<5+27 



1 1) (1 1 . The effect of the previous LC res- 



onator can be neglected because it is far off resonance with 
the PCB. However for the purpose of universal QIP this sec- 
ond LC resonator is not absolutely necessary |i2g| , as will be 
shown below. 

Inductively coupled PCB arrays. - To realize universal QIP 
on the PCBs, we need a scheme to couple them. We notice 
that different PCBs will experience different noise as they are 
biased by different leads. Those far apart are susceptible to 
different charge noise too. So we only have what we call "lo- 
cal" DFS; this only relies on noise symmetry over a few, here 
two, physical qubits, as is inevitably the only realistic case for 
scalable QIP, and previously discussed methods |20] do not 
apply. Stringent restrictions are put on the two qubit coupling 
in order to preserve the DFS. A capacitive coupling between 
lb and 2a, for example, cannot be used, as this would cause 
the noise in PCBl's lead to leak asymmetrically into PCB2. 
Furthermore neither PCB may leave its DFS during its evolu- 
tion. 

Here we discuss a new approach that allows scalable QIP 
based on local DFS. We couple the PCBs inductively, using 
a small mutual inductance AI' between their dc-SQUIDS, as 
shown in Fig. [T](b). As the dc-SQUIDS are biased at $o/2 
the coupling Hamiltonian is AcrJ^jCrf^ l23l . where the coupling 
strength A = j^I' Ip. ib^c 2a (chosen much smaller than E'j, 
the unsuppressed Josephson Energy of the SQUIDS). 

Obviously the PCBs in the array in Fig. ^(b) can be ini- 
tialized in their DFSs just as described before. Now, if all the 
PCBs are biased at the same voltage, they will get out of their 
DFSs due to resonant interactions by the above coupling term. 
Therefore, we bias the odd numbered PCBs at the degenerate 
point Ug ^ 1/2, giving a = and the even numbered 
PCBs at some other point (but far off resonance with their LC 
resonators) giving a large i?^. Because of the large detuning, 
B^ 3> A, the coupling Hamiltonian has no effect in the idle 
mode lll2il . One might worry that this different biasing will 
introduce a difference between the phase frequencies of the 
PCBs, but it is not the case since the DFS states always have 
the energy 7 regardless of the voltage bias of the PCBs. When 
we want to do a controlled phase gate between PCB 1 and 2, 
we tune their biases near some common target value very dif- 
ferent from their previous values (such that they do not inter- 
act with other neighboring PCBs) and the LC frequencies. As- 
suming the fields are Bf and we work in the dispersive re- 
gion such that the states of the PCBs do not change except that 
dispersive phases are obtained due to the virtual transitions. 
For instance, the energy of the state IO1O2) will be shifted by 

^"^:^;^,!:;'^ - where \m) ^ | UMbhM is 

the virtually excited intermediate state and A = Bf — B2 
is the detuning. Other phases can be calculated too, giv- 
ing in the basis IO1O2), IO1I2), II1O2), II1I2) an effective 
Hamiltonian diag{j^, 47+(B?+i3|) ' 47-(4'+i3|) ' 4^1' 
which reduces to diag{Q^ 0, 0, } in the dispersive region 
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A <C |47-A| < |47 + A| < |47± (Sf + This gives a 
CPHASE(a) gate diag{l, 1,1, e~*"}, where a = with 
t being the evolution time. 

Notice that in the above implementation of the CPHASE 
gate, the capacitive coupling 7 plays an essential role. As can 
be easily checked the above procedure does not give an en- 
tangling gate if a and b do not bias each other (7 = 0). This 
is because IO1O2} and II1I2} would acquire opposite energy 
shifts ±A^/A, due to the fact that the energy differences be- 
tween the initial and intermediate states are opposite for these 
two cases. Thus the role of the capacitive coupling 7 can be 
understood by an interesting "parity argument;" under the ex- 
change of the states of a and b for both PCBs, the energy mis- 
match (the denominator in the perturbative calculation) due to 
the detuning A is odd. This symmetry is broken by the cou- 
pling between a and b, which is even under this operation (as 
is obvious from the form — jcr^a^). 

Another point of interest is that, a phase gate on a PCB 
can be implemented by using the CPHASE(q:) and the X gate, 
as is easily recognized by the identity e'^^" = e^"{X2 ■ 
C P H AS E {2a))^ . Therefore, for the purpose of universal 
QIP on the PCB array the LC resonator giving the phase gate 
is not absol utely necessary, as long as the system has more 
than 1 qubit I2al . This is potentially beneficial in reducing the 
hardware. 

Discussion. - In the above we described our scheme of QIP 
with PCBs based on encoded qubits and dispersive dynamics, 
which requires only tuning the voltage bias of the PCBs. Our 
scheme can prevent decoherence from collective noise. The 
circuit noise obviously couples symmetrically to the CBs of 
the PCB. The charge noise requires some caution, since its 
exact nature is still a topic of debate fl7ll27ll . The early ex- 
periments in flTil clearly show that the charge noise on close 
by islands are correlated. The conclusion drawn from this ob- 
servation, that the charge noise stems mostly from sources in 
the substrate was further substantiated by 12811 . This has im- 
portant consequences, because it suggests that it is possible 
to engineer the environment for desired noise configurations. 
Indeed, analysis in 1 17] shows that high noise correlations can 
be achieved for properly designed geometry and layouts of the 
charge islands. Simple environment engineering was already 
successful lEslEgil in various contexts. 

For our scheme to work, the CBs must be located within 
a distance smaller than the wavelength of the background 
charge fluctuations, so that they experience the same noise. 
This seems to be realistic, since the advance in device fab- 
rication allows to make smaller structures, and more impor- 
tantly the fact that the charge noise originates from the sub- 
strate makes it possible to engineer the environment for the 
desired noise symmetry ITtIEsIi . For instance^if we put the 
PCB on an electrode instead of the substrate 12811 . the charge 
impurities will be located far away from the PCB, and thus 
couple symmetrically to the CBs. 

Though our scheme eliminates the effect of the collective 
charge noise on the PCB, the decoherence time will be finite 
as there are other non-collective noise in the system not dealt 
with by our prescription, for instance dissipation due to the 
finite impedance of the junction and the noise in its critical 



current. The effect of the virtually excited states and the fluc- 
tuation of the dispersive energies must be evaluated carefully 
too, though some results were obtained previously l30l I3TI1 . 
Qualitatively, as shown by simple analysis based on Master 
equations the number of operations allowed in our dispersive 
scheme increases hy A/g ^ 1 (here, A the effective detuning 
and g the coupling strength) as compared to the usual scheme 
based on resonant Rabi manipulations, if the same coupling 
strength g is assumed 1 32 j. The PCBs do not experience de- 
coherence in the idle mode, which is favorable for a large sys- 
tem in which only a fraction of the qubits undergo active ma- 
nipulation at the same time. Therefore, our scheme we can 
reduce the error rate of the PCBs below the threshold for error 
correction schemes l33ll and thus make superconducting QIP 
feasible. Detailed calculation of the decoherence time for a 
realistic PCB system will be reported elsewhere i32il . 

Another practical concern is that the CBs in a PCB will not 
be completely identical due to the imperfect fabrication. Be- 
cause only local symmetry is required, this problem is less 
significant since fabrication variations tend to happen at large 
scales and experiments show that closely spaced charge-boxes 
can be equal beyond experimental resolution Il7i1 . In addi- 
tion, the error induced by non-identical qubits was shown to 
be higher order in the symmetry breaking ■ Therefore, we 
conclude that the technological problem of imperfect fabrica- 
tion is already solved to the extent needed for our scheme. 

Parameters. - Finally, we give some parameters for the ex- 
perimental consideration. We use small CBs closely spaced 
with a total capacitance Ct « 0.16/F and charging energy 
Ec ~ 500GHz. A mutual capacitance Cc = 5aF gives 
7 « 7.5GHz. A mutual inductance M = 7pH between the 
PCB and LC with L = 50pH and uj/2n = 200GHz gives 
g = Q.25GHZ for 1^ = iOnA. Tuning the bias of the PCB 

close to the LC frequency with a detuning 5 « —12.5GHz 

2 

results in an exchange interaction with the strength ~ 
2hMHz, corresponding to a period of 40ns. During the idle 
mode, we bias the odd numbered PCBs at = and the 
even numbered ones at B~ k, lOOGHz. Since low self in- 
ductance dc-SQUIDS should be used, we choose M'lc ~ 
10~^(l>o nil' and M' ^ 120pH gives the coupling strength 
A w 0.22GHz. Tuning the biases of the neighboring PCBs 
both to about AOOGHz with a detuning A « 28GHz gives a 
CPHASE gate at the rate w 25MHz. The above pa- 

rameters are well within the reach of the current technology 

El. 

In conclusion, we have discussed a technique for robust and 
easy superconducting QIP. By combining the ideas of encod- 
ing and dispersive manipulations, we protect the charge qubits 
from the dominating decoherence and realize universal QIP 
on the encoded qubits with minimal control. Besides the great 
potential of solving the fundamental difficulties in supercon- 
ducting QIP, we expect the general idea of dispersive manip- 
ulation of encoded qubits to be of interest to other physical 
systems, such as atomic and other solid state systems ISol . 
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